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Abstract 

Gauge theory is a well-established concept in quantum physics, 
electrodynamics, and cosmology. This theory has recently proliferated 
into new areas, such as mechanics and astrodynamics. In this paper, 
we discuss a few applications of gauge theory in finite-dimensional 
dynamical systems with implications to numerical integration of dif- 
ferential equations. We distinguish between rescriptive and descriptive 
gauge symmetry. Rescriptive gauge symmetry is, in essence, re-scaling 
of the independent variable, while descriptive gauge symmetry is a 
Yang-Mills-like transformation of the velocity vector field, adapted to 
finite-dimensional systems. We show that a simple gauge transfor- 
mation of multiple harmonic oscillators driven by chaotic processes 
can render an apparently "disordered" flow into a regular dynamical 
process, and that there exists a remarkable connection between gauge 
transformations and reduction theory of ordinary differential equa- 
tions. Throughout the discussion, we demonstrate the main ideas 
by considering examples from diverse engineering and scientific fields, 
including quantum mechanics, chemistry, rigid-body dynamics and 
information theory. 



1 Introduction 



In modern physics, gauge theories are probably among the most pow- 
erful methods for understanding interactions among fields. The im- 
portance of gauge theories for physics stems from the tremendous suc- 
cess of the mathematical formalism in providing a unified framework 
to describe the quantum field theories of electromagnetism, the weak 
force and the strong force. Modern theories like string theory, as well 
as some formulations of general relativity, are, in some sense, gauge 
theories. The Yang-Mills theory, the standard approach to quantum 
field theory, is a particular example of gauge theories with non-Abelian 
symmetry groups. Gauge symmetries are the core mathematical mech- 
anism of gauge theory, reflecting a redundancy in a description of a 
system. 

The earliest physical theory which had a gauge symmetry was 
Maxwell's electrodynamics. However, the importance of this symme- 
try remained unnoticed in the earliest formulations. After Einstein's 
development of general relativity, Hermann Weyl, in an attempt to 
unify general relativity and electromagnetism , conjectured that "Eich- 
invarianz" or invariance under the change of scale (or "gauge" ) might 
also be a local symmetry of the theory of general relativity. Weyl 
coined the use of gauge symmetry in modern physics j4Ul 1411 H^ . say- 
ing that "symmetry, as wide or narrow as you may define its meaning, 
is one idea by which man through the ages has tried to comprehend 
and create order, beauty, and perfection". 

The question to be raised at this point is: Can the gauge-theoretical 
approach be applied on finite-dimensional systems, and to what de- 
gree of success? Keeping in mind that Lie symmetry has been a major 
tool in the study of ordinary differential equations (ODEs) [21], the 
answer must be positive; however, we have not established yet a clear- 
cut connection between Lie symmetry and gauge symmetry. We plan 
do so in the sequel, and show that, indeed. Lie point-symmetry is 
closely connected to gauge symmetry, albeit this connection is not al- 
ways straightforward. In order to facilitate the establishment of such 
a connection, we use gauge symmetry in a more general context, a 
context of a symmetry defined by diffeomorphisms. This will ulti- 
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mately allow us to combine various manifestations of gauge symmetry 
in finite- dimensional dynamical systems under a single mathematical 
realm. 

The idea to apply gauge transformations on ODEs is not new; it 
was suggested by Kunin j2L)j . and has been recently revived by Efroim- 
sky El El IS21 E3 El ) who developed a gauge-generalized astro- 
dynamical theory for modeling the effect of orbital perturbations using 
non-osculating orbital elements. Eforimsky's gauge theory, however, 
does not deal with scale transformations, while Kunin's gauge theory 
is mostly concerned with local gauge groups and discrete symmetries. 
Thus far, there has not been a unified gauge theory which is able to 
support both Yang-Mills-like gauge theories and scaling theories in 
finite-dimensional systems. 

In the current paper, we shall attempt to fill this gap by developing 
a gauge theory for finite-dimensional dynamical systems through two 
gauge symmetry mechanisms; the first symmetry mechanism will be 
called rescriptive gauge symmetry, evoked by carrying out a rescriptive 
gauge transformation. 

Rescriptive gauge symmetry succumbs to the fundamental notion 
of gauge transformations, namely, a change of scale, and is also inti- 
mately connected to bilinearity. To show rescriptive gauge symmetry, 
we shall carry out an infinitesimal transformation of the independent 
variable - which in the bulk of our subsequent discussion will be time 
- into a different scale. The manifestation of gauge symmetry in this 
case will be reflected in the ability to obtain equivalence between the 
direction fields of the original and gauge-transformed systems. In 
many practical applications, this implies that the system can be re- 
duced to a form amenable for quadrature (e. g. linear ODEs). We shall 
formalize this observation by establishing a mechanism for reduction 
through rescriptive gauge symmetry. 

To illustrate the concept of rescriptive gauge symmetry, we will 
present a myriad of physical examples taken from diverse scientific and 
engineering fields, including rigid-body dynamics, finite-dimensional 
quantum mechanical systems, chemistry, and information theory. 

An instrumental constituent of our new theory is the gauged pen- 
dulum. Generally speaking, a gauged pendulum is a physical system 
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with a quadratic integral of motion, whose behavior in the time do- 
main can be arbitrary, although its phase space structure remains in- 
variant under a change of scale. This implies that after a suitable scale 
transformation, harmonic oscillations will emerge. We show that many 
physical systems can be either re-formulated to match the formalism 
of the gauged pendulum, or are natural gauge pendulums per-se; a 
classical example for a natural gauged pendulum is the Euler-Poinsot 
system, to be subsequently analyzed. 

We ultimately utilize the notion of a gauged pendulum to question 
some common engineering misconceptions of chaotic and stochastic 
phenomena, and show that seemingly "disordered" (deterministic) or 
"random" (stochastic) behaviors can be "ordered" , or, put differently, 
evoke simple patterns j281 129j using an infinitesimal transformation of 
the time scale. This brings into play the notion of observation and 
observables; we show that temporal observations may be misleading 
when used for chaos detection. 

The second symmetry mechanism, reminiscent of the gauge sym- 
metry arising in Maxwell's equations [HJ ^-nd its generalization into 
the Yang-Mills field theory, will be referred to as descriptive gauge 
symmetry. The concomitant gauge transformation will be called a 
descriptive gauge transformation. Descriptive gauge symmetry natu- 
rally arises in Newtonian mechanical systems, and can be thought of as 
an invariance of some configuration space under a gauge transforma- 
tion of the covariant derivative. In fact, descriptive gauge symmetry 
may be best understood by relating it to the method of variations- 
of-parameters (VOP), which is an analytical formalism for solving 
inhomogeneous (forced) differential equations. 

Euler invented the VOP method jl4l I15j for treating highly non- 
linear problems emerging in celestial mechanics. However, it was La- 
grange who employed this method for deriving his system of equations 
describing the evolution of the orbital elements |211 1221 123j . known 
as Lagrange's Planetary Equations. The relation between the VOP 
method and descriptive gauge symmetry can be explained as follows. 
According to the VOP method, the integration constants of the homo- 
geneous solution of a given ODE are endowed with a time variation 
due to the presence of an external force. However, the transforma- 
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tion from the state variables of the original problem's phase space to 
the new state variables defined as the time- varying constants involves 
an inherent freedom, which, in practical calculations, can be removed 
by means of a user-defined constraint. The constraint may be essen- 
tially arbitrary insofar as it does not come into contradiction with 
the equations of motion written for the variable "constants." The in- 
ternal freedom emerges under the following circumstances: First, one 
should perturb some n-dimensional differential equation, and solve 
it by the VOP method (i.e., using the unperturbed generic solution 
q(t, xi . . . Xn) as an ansatz, and making its constants Xj time- varying) ; 
second, the number of "constants" promoted to variables must exceed 
n. Thus, when the said equations are written as equations for the new 
state variables Xj, the number of these variables will exceed that of 
equations; hence the internal freedom. Mathematically, this freedom 
is analogous to the gauge symmetry in electrodynamics, while the re- 
moval of this freedom by imposing an arbitrary constraint is analogous 
to fixing of a gauge in the Maxwell theory. 

Prom a practical standpoint, we show that descriptive gauge sym- 
metry may be used to considerably mitigate the numerical truncation 
error of numerical integrations, and even "symplectify" non-symplectic 
integrators. Wc also show how a given system of ODEs can undergo 
a reduction under a descriptive gauge symmetry transformation. 



2 Preliminaries and Definitions 

Consider a finite-dimensional dynamical system whose dynamics are 
modeled using first-order vector differential equations of the form 

where x G M", f : M" R". Assume that this system is subjected 
to structural modifications resulting from some re-formulation of the 
phase space, perturbations, control inputs, exogenous disturbances 
or modeling uncertainties. We shall generalize these modifications 
under a single mathematical umbrella which we call a rescription. A 
rescription operator in the time domain, ^t, acts upon the vector field 
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f in the following manner: 

S^t of = g(x,u(x, t)) 



(2) 



The vector field u : x M ^ MJ^, m <n will be called a rescriptor. 
The rescriptor may be either static or dynamic. In the former case, 
one may write u = K{x, t), with K : M" x M ^ M™, while in the latter 
case the rescriptor constitutes a dynamical system of the general form 



where y G C M^, M. being some compact differentiable manifold, 
and hi : X — is a vector field. By restricting the dynamic 
rescriptor from dependance on u, we can unify the dynamic and static 
rescription, viz. a static rescription can be defined as a special case of 
a dynamic rescription. 

In most cases, the rescription modifies - sometimes intentionally, 
such as in the case of control inputs - the fundamental properties of 
the original system. These "fundamental properties" may be, for in- 
stance, integrability, symmetry and structure/ volume-preserving mea- 
sures. The new properties of the rescribed system can be investigated 
both in the time domain and in the phase space. 

However, in some cases the rescription is merely an illusion; that is, 
the rescription does not change the phase space and the fundamental 
properties of the original system, although it could modify the flow 
ip{xi{t = t()),t). A classical example is the action of the rotation group 
Go = SO{n) on Lagrangians of the form C = /cx^, which remain left- 
invariant under the transformation x ^ Tx, T € SO{n). We shall 
exclude these trivial occurrences from our discussion, and will explore 
a more general setting. In this general setting, the system is invariant 
under the action of some (possibly time-varying) finite-dimensional 
gauge group, Q. 



du 

dy 
dt 



hi(y,u). 



h(x, (ix/(if,y, u) 



(3b) 



(3a) 
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3 Rescriptive Gauge Symmetry 



We ask whether the system can be "de-rescribed" by finding new inde- 
pendent variables, Tj, possibly different for each rescriptor component 
Ui, satisfying 

dTj = Gi{x, Ui(x, t), dx, dt) (4) 

for which 

5, o g = f (x) = x', (5) 

where the operator ()' denotes differentiation of each xi with respect 
to some Tj, 

x' = ^, i = l,...,n, je[l,...,n]. (6) 

dTj 

If 3(iTj, j G [1, . . . , n] satisfying such that © holds, then we shall 
say that system (j2| exhibits full rescriptive gauge symmetry^ under 
the rescriptive gauge transformation @. In this case u becomes either 
a static or a dynamic rescriptive gauge function. 

A rescriptive gauge symmetry of order p or simply partial rescrip- 
tive gauge symmetry comes about when the rescriptive gauge trans- 
formation de-rescribes only p state variables, p < n, viz. 

drOgi = f^{x) = Xi', iGW. (7) 

In this case, if t S t = [0, tf),tf < oo and Tj G R, j G [1, . . . then 
3Tjo, to,Xi(to),Xi{TjQ) such that the flow satisfies the gauge homeomor- 
phism 

(p{Xi{to),t) = (p{Xi{Tjo),Tj) (8) 

for t n R, where the flow is interpreted as the one-parameter group of 
transformations 

Gt : Xi{to) ^ Xi(t), Gr ■■ Xi{To) ^ Xi{T). (9) 

The notion of rescriptive gauge symmetry has far-reaching applica- 

^The definition of rescriptive is given in Webster's Revised Unabridged Dictionary 
(1913): "Pertaining to, or answering the purpose of, a rescript; hence, deciding; settling; 
determining" . 
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tions in "ordering" seemingly "disordered" phenomena, solving ordi- 
nary differential equations (ODEs) and improving numerical integra- 
tion thereof. We shall illustrate these ideas by discussing a few exam- 
ples of practical interest. We embark on our quest by presenting the 
notion of a gauged pendulum, dwelt upon in the following subsection. 

3.1 The Gauged Pendulum 

Finite-dimensional systems can often be modeled by Hamiltonian vec- 
tor fields induced by a nominal Hamiltonian, Ti., and a perturbing 
Hamiltonian ATC. Moreover, in ubiquitous fields of science and en- 
gineering, TC is comprised of n uncoupled harmonic oscillators [2], 
namely 



1=1 

n 

= ^T.n.Mt),Mt)] (10) 

i=l 

where Q, = diag(a;f , . . . , w^), 

c[=[qi,...,qnf, p=[pi,...,pnf, (11) 

are the generalized coordinates and conjugate momenta, respectively, 
so that (q, p) G T*Q, where Q is the configuration space, T*Q is the 
cotangent bundle of Q, and dim T*Q = 2n is the dimension of the 
phase space. 

Hamilton's equations for i = 1, . . . ,n are then 

Qi = Pi (12a) 
Pi = -Lofqi. (12b) 

Carrying out the point transformation into action-angle variables, 
given by 

Qi = \ — sin (pi, Pi = \f^^j cos (13) 
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simplifies the Hamiltonian ()1U|) even further, into 

H[u:{t),m] = = ^f;u;,<I>„ (14) 

1=1 

where 

uj = [u;u...,u;nf, * = [$1,...,$,]^. (15) 

We note that from the topological standpoint, in both (|10jl and (|14)1 
7^ G S*" C M"^^ is always homeomorphic to an n-ellipsoid, and Tli G 

is an integral. Moreover, the transformation (q, p) i-^ {(p, $) is a 
(universal) covering map of S"'. 

It is clear that the Hamiltonian (|10|) is left-invariant under the 
action of the rotation group Go = SO{n) on q and p if = wq- 
However, we shall seek a broader invariance of 7i with respect to 
the gauge group, which does not necessarily adhere to the SO{n) 
symmetry. 

To that end, let us choose an arbitrary (not necessarily smooth) 
scalar field tij(q, p) : T*Q — > M to serve as our rescriptor, coupling 
the dynamics of the n pendulums, and re-write (|12)) into the strictly 
bilinear form^ in [q,p\ and u: 

Qi = PiUi{q, p) (16a) 
Pi = -uJiQiUiiq, p). (16b) 

Obviously, a constant of motion for each of the pairs {qi, pi) would be 

Ci = ^{p^+u;fqf), i = l,...,n (17) 

although Ci is no longer the Hamiltonian. Nevertheless, system H16|l 
remains integrable regardless of the particular form of Uj, since there 
are n integrals for n degrees-of-freedom. This can readily observed by 
performing the (affine in dt) rescriptive gauge transformation 

dn = Ui{p,q)dt, (18) 

strictly bilinear system with respect to x and u has the structure x = Mxu, M g 
R"^". See for details. 
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which, on one hand, extends ()16|) into the state-space model 



qi = PjU.j(p,q) (19a) 
Pi = -iofqiUi{i>,(i) (19b) 
fi = Ui(p,q), (19c) 



but, on the other hand, transforms (|16|) back into the simple har- 
monic oscillator form in the independent variables Tj, assuming the 
symplectic structure 

q'i = Pi (20a) 
Pi = -^hi- (20b) 

Thus, Ui is a rescriptive gauge function Vi, and Cj can be interpreted 
as the Hamiltonian again, that is, 

W[q(r,), p(r,)] = ^ = I^Y^^pI + -^hh ■ (21) 

i=l i=l 

In this example the transformation t i— > rj,i = 1, . . . ,n is therefore 
a static rescriptive gauge transformation. This means that Ui may 
be used to control the flow of pi and qi in the time domain, but the 
persistence of the integrability under the transformation H18|) forces 
the system to exhibit the same behavior as the harmonic oscillator in 
the modified times for each degree-of-freedom. 

In order to generalize this concept and illustrate how rescriptive 
gauge functions emerge in common physical systems, we must allow 
Ui to be an output of a dynamical system, giving rise to dynamic 
rescription, defined in ^ In this case, system (|lf)j) . written for each 
degree-of-freedom, z = 1, . . . , n, becomes 

qi = PiMi(p,q) (22a) 

Pi = -ujfqiUi{p,q) (22b) 

iii = /ii(p,q,Ui,y) (22c) 

y = hi(u,y). (22d) 

Carrying out the rescriptive gauge transformation H18|l reveals a par- 
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tial rescriptive gauge symmetry: 



Qi = Pi 

Pi = -^Ui 

u'i = — /ii(p,q,Mj,y 



(23a) 
(23b) 
(23c) 



Thus, independently of the particular characteristics of the dynamic 
(or static) rescriptive gauge function, Ui, the system re-assumes the 
harmonic oscillator structure for {qi,Pi). This situation can therefore 
be viewed as a generalization of the pendulum model. The persistence 
of the harmonic oscillations under the rescriptive gauge transforma- 
tion gives rise to the concept of a gauged pendulum. The gauged 
pendulum is a dynamical system whose flow becomes periodic under 
the rescriptive gauge transformation, although the flow of the origi- 
nal system may exhibit arbitrary behavior in the time domain. Such 
systems arise in ubiquitous fields of science and engineering. For ex- 
ample, the following model arises in the study of quantum mechanical 
phenomenon (assuming a zero decoherence coefficient) |18) : 



This is obviously a gauged pendulum with the static rescriptor u = 
uiU2- we shall subsequently dwell upon additional physical examples. 

An alternative formulation of systems exhibiting partial rescrip- 
tive gauge symmetry with a dynamic rescriptive gauge function may 
written as 




(24) 
(25) 







(26b) 
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which becomes 



di 
P'i 



Ui = 



Pi 

-^hi 

— /ii(p,q, Ui 

Ui 



(27a) 
(27b) 

(27c) 



after de-rescription using our standard rescriptive gauge transforma- 
tion. Here the rescriptor Ui still constitutes a dynamic rescriptive 
gauge, albeit it is not an output of an auxiliary dynamical system 
anymore. In fact, if we relieve hi from direct dependance upon Uj, viz. 
Ui = hi{-p, q), then we uncover additional integrals of the motion, /Cj, 
defined by the quadrature 



K-i 



1 



:Ui 



J hi[ci{Ti),p{Ti)]dTi. 



(28) 



These new constants posses a clear meaning, revealed by writing 

die. 



Ui 



On 

dlCi _ 
dui 



hi 



Ui 



(29a) 
(29b) 



Hence, Ti,Ui can be interpreted as generalized coordinates and conju- 
gate momenta, respectively, evolving on a 2n-dimensional symplectic 
manifold. /Cj is then a Hamiltonian, and the dynamics of (uj,rj) is 
integrablc. This remarkable structure implies that if the time deriva- 
tive of the rescriptor docs not explicitly depend upon the rescriptor, 
then the rescriptive gauge transformation may be viewed as the sym- 
plectomorphism 



_ dniqi{t),pi(t)) 
d^i 

dn{qi{t),piit)) 

P'- Wi 



Pi 



dH{qi{Ti),pi{Ti)) 

_dTi{qi{Ti),pi{Ti)) 
^i 



Ui 



dn 



dJCi 
dui 



(30) 
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The formulation in ()3U|) is general, and is not limited to Hamiltonians 
of the form (|21|1 : rather, if Ti.i{pi, qi) = const, is a given Hamiltonian, 
then Hamilton's equations 

. _ dni{pi,qi) . _ dHiipuqi) , . 

' Fq ^^^^ 

undergoing a rescription 

dHiipi^qt) , . . dHi{pi,qi) , . . 
Qi = ^^i(P,q), Pi = Ui{P^W (32) 

will still possess still TCi{pi,qi) = const, as an integral, and can be de- 
rescribed using the rescriptive gauge transformation dr = Ui{p,q)dt 
into 

, dni{pi,qi) dni{pi,qi) , , 

'^^ = dp ' P= 

Finally, a slightly different formulation of the gauged pendulum with 
a dynamics rescriptive gauge symmetry, to be illustrated in ^3.3| may 
be written as 

(ji = kiPiUi{p,q) (34a) 
Pi = k2qiUi{p,q) (34b) 



Ui = hi{p,q,Ui) (34c) 



which becomes 



q'i = hiPi (35a) 
Pi = hiqi (35b) 
u'i = —hi{p,q,Ui) (35c) 

Ui 

after de-rescription using the rescriptive gauge transformation = 
Uidt. Here we have 

= ^ {hipf - k2iqf) , i = l,...,n (36) 

as integrals. However, an important caveat is that (|35a)) - ()35b|) may 
be viewed as a gauged pendulum only if kiik2i < 0. Otherwise, the 
de-rescription will yield hyperbolic motion in the variable r. 
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3.2 Newtonian Systems 

We shall now show that Newtonian systems can be re- written into the 
gauged pendulum formalism. To that end, consider the system 

q = vp (37a) 
p = —vq (37b) 

V = f{x,v). (37c) 

Performing the "action-angle" transformation (cf. Eq. (|13|)) p = 
cosx, q = sin a; and re- writing 1)3 7|) yields 

X = V (38a) 

V = f{x,v), (38b) 

which is a state space representation of the Newtonian system 

x = f{x,x). (39) 

We immediately observe that our rescriptor is the velocity, v. Hence, 
when dealing with problems in the Newtonian context, the rescriptive 
gauge function is simply the gauge velocity. 

Our Newtonian system therefore possesses a trivial partial rescrip- 
tive gauge symmetry, found by performing the rescriptive gauge trans- 
formation dr = dx = vdt. In other words, in the Newtonian case 
r = X, and the de-rescribed system becomes 

q' = p (40a) 
p' = -q (40b) 
v' = f{x,v)/v, (40c) 

where ()' denotes differentiation with respect to x. If f{x,v) = f{x), 
then ^ 

^ = \- j f(^)dx (41) 

is an integral, K. is the Hamiltonian for the original system (|38j) . and 
Tl = {q^ +p^)/2 is the Hamiltonian of {q,p)- Thus, any Newtonian 
system can be written in the gauged pendulum form by extending 
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the phase space dimension by one. Consequently, any system whose 
state-space model is similar to (jHTj) is a Newtonian system in disguise. 



3.3 Eulerian Systems 

In a body-fixed frame, the attitude dynamics of a rigid body are usu- 
ally formulated by means of the Euler-Poinsot equations. In a free- 
spin case, these equations look as 

lu + u X lu = 0, (42) 

I being the inertia tensor and u = [wi, W2,W3]"^ G © is the body an- 
gular velocity vector, where & is the foliation {{luji, Iuj2, /tJ3)|IiLi;^ + 
l2t^2 + Is'-^s = G*^}), G being the total angular momentum. 

Assuming that the body axes coincide with the principal axes of 
inertia, 

I = diag(Ii,l2,l3), (43) 
the Euler-Poinsot equations are 

uji = aiuj2^3 (44a) 
i02 = (T2iOii03 (44b) 
= cj3u;iu;2 (44c) 

where 

h — h h — h h — h 

0-1 = r , 0-2 = r ' ^5 = F • (4^) 

h h ^3 

We shall now show that the Euler-Poinsot equations are a classical 
example of the gauged pendulum concept with a dynamic rescriptive 
gauge, exhibiting partial rescriptive gauge symmetry of order 2. To 
that end, define the rescriptive gauge transformation 

dr = ojsdt (46) 
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and re-write ()44() into 



uj'i = a\uJ2 (47a) 
= 02^\ (47b) 
= — a;iu;2, (47c) 

which adheres to the gauged pendulum model H35|) . Thus, in the 
modified scale r, wi and will exhibit harmonic oscillations with 
frequency yjfJio^ if crio-2 < 0, given by 

-0-2^10 sin(a;oT"o) + ^20^^o cos(a;oTo) , . 
W2(t) = cos(l<jo''") 

Wo 

cr2a;iocos(woTo) +^^20'^osin(6Jo-ro) . , . 
+ sm(a;oT) (48) 

uo 

(T2a;io cos(woTo) + t*^20'^o sin(cJoTo) , , 
ti;i(r) = cos(wo''") 

W2oWocos(a;oT"o) - 0-2^10 sin(a;oT-o) . , . , 
— sm(iwo''") (49) 

^2 



where = vkio^, tJio = wi(ro), a;2(ro) = t^20- 

The solution for the dynamic rescriptor can now be easily solved 
by quadrature.^ Since 

C = ^t^i ~ ^3 y t^ii:J2C?T (50) 

is an integral, 

= -\/2C + A(T3 cos2(a'o''") + i?o"3 sin(a'oT) cos(a;oT) (51) 



3 



Note that here the rescriptor has units of angular velocity, while in the Newtonian 
case it was the velocity. We shall re-iterate on this issue in the following sections. 
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where 
A 



-2cj2 cos^(u;o''"o)i^io ~ 4cos(a;oTo)a;io'^20 s\ii{ujQTQ)uJoa2 



+ 



2, ,2 
10 



'^0'^2 



(52) 



B 



-2a2 cos{ujoTq)ujIq sin(a;o''"o) + 4cos(a;oTo)^cdiot<J20i^of2 



"^0^2 



+ 



2a;2o'^o'72Wio + 2^0^20 sin(wo''"o) cos(a;oro) 

^0^2 



(53) 



From (|4H)1 . the new independent variable is 



(54) 



To understand its physical meaning, we recall that if we let the hat 



map 



so (3) denote the usual Lie algebra isomorphism that 



identifies (5o(3), [, ]) with (M^, x), then 



where 





UJ = 


-RR^ 













UJ = 







-LJl 




—UJ2 








(55) 



(56) 



and R G SO (3) is a rotation matrix from inertial to body coordinates. 
If we take the rotation sequence 3— >1— >3, c/)^^— evaluation 
of (|55|) will entail the well-known expressions for the components of 
the vector of the body angular velocity u in terms of the Euler angles 
rates (/>, 9, ip: 



voi = (f) sin 9 sin x/j + 9 cos ip 
UJ2 = (j) cos ip sin 9 — 9 sin tp 
= ■0 + cos 9 



Thus, 



UJ3 = 1p + 



> COS 9dt, 



(57a) 
(57b) 
(57c) 

(58) 



17 



so the coi — T dynamics may be viewed as a solution for the phase space 
of the Eulerian system. 

The rescription in the Eulerian case possesses an interesting sym- 
metry. In the above discussion, we detected the dynamic rescriptor 
^3 for the pair (wi, L02), but there are two other possible rescriptions: 
LOi for {uj2, UJ3) and uj2 for ("^ii ^3)- 

Finally, we note that the rescriptive gauge transformation lin- 
earized the Euler-Poinsot equations; observe that H47a() - H47b|l are lin- 
ear, and H47c() is a simple linear quadrature thereof in the variable 
z = w|. We shall further dwell upon this finding in ^l-j.HI 



3.4 Other Common Systems Exhibiting Re- 
scriptive Gauge Symmetry 

The gauged pendulum is a particular case of systems exhibiting re- 
scriptive gauge symmetry. However, there are systems that exhibit 
rescriptive gauge symmetry, which cannot be rendered periodic after 
a rescriptive gauge transformation. Generally speaking, such systems 
cannot be conveniently described using the Hamiltonian formalism, 
although they do posses integrals. Consider, for illustration, the dy- 
namical equations of two chemical reactants, A and B, whose concen- 
trations evolve according to the bilinear rate law [5] 

"^^^^ - ki[A][B] (59a) 



dt 
d[B 

~dt 



k2[A][B]. (59b) 



These can be de-rescribed using e. g. dr = [A\dt, yielding the linear 
equations 



dr 

d[B] 

dr 

so that 



ki [B] (60a) 
k2[B], (60b) 



[Bir)] = [Biro)]e'-\ [^(^)] = [S(ro)]|i(e'=- - 1) + [^(to)]. (61) 

k2 
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An integral for system ()6U() is C = [A] — ki/k2[B] , albeit this is not the 
Hamiltonian. Consequently, an additional class of systems exhibiting 
rescriptive gauge symmetry may be written as 

Qi = qiUi{p,q) (62a) 

Pi = qiUi{p,q) (62b) 

iii = /ii(p,q, Ui,y) (62c) 

y = hi(u,y). (62d) 



3.5 The One-Parameter Lie Symmetry Group 

Thus far we have not explicitly spelled out a relationship between the 
rescriptive gauge transformation and Lie point-symmetry transforma- 
tions. This is the purpose of the following discussion. 

To keep things simple, assume a 1-DOF gauged pendulum model 
with a static rescriptor, u{p, q): 

q = pu{p,q) (63a) 
p = -qu{p,q). (63b) 

This set of equations can be analyzed by means of one-parameter 
groups based upon infinitesimal transformations. We demand the 
equation to be invariant under infinitesimal changes of the indepen- 
dent variable t, but without a simultaneous infinitesimal changes of 
the dependent variables. This leads to the Lie point-symmetry trans- 
formation 

p ^ p (64a) 
q ^ q (64b) 
t ^ T = t + eC{p,q). (64c) 

We now apply (|6l|) on (|63|) by following these stages: First, we write 

^'^ +0(.2) (65) 



dt + e['it,dq+'t,dp 
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and 

T- = ^ ^ + 0(^'). (66) 

Expanding and (|^^ into a Taylor series with e as a first-order 
small parameter we get 

dq dq dq / d( dq dC dp\ ^ 
d^ = Tt-'Yt[d-qTt^d-p^)^^^'^^ 



dp dp dp f dCdq dC dp\ 2n /«o^ 

d^=d^-'^[d-qd-t^d^d^)^^^'^- 

These yield a partial differential equation (PDE) for C{p,q), 

u{p,q)e\p — — q — — — j - u{p,q) + 1 = 0, (69) 

the solution thereof is 

<(P, ,) = - r ^WEZ>=i=^., + c„c (70) 

where c = and cq is an integration constant. To relate (|65|) and 

()66|1 to the generators of the infinitesimal transformation, we write 

T = t + eC{p, q) + ... = t + eXt + ... (71) 

where the operator X is given by 

X = C(P,9)| (72) 

In addition, due to the fact that ()63() is autonomous, it will also exhibit 
Lie point-symmetry with generator 

X. = |. (73) 

Symmetries (|72|) and (|7,3|1 form an Abelian Lie algebra X with the Lie 
bracket [X, Xi] = 0. 
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In essence, this symmetry implies that the direction field is 



dp 
dq 



P 



q 



(74) 



and is therefore homogenous, that is, invariant under all dilations 
{p, q) ^ {e^p, e^q), A G M, which holds true for any dynamic or static 
rescriptive gauge u{p, q). The connection to the rescriptive gauge sym- 
metry can now be easily obtained via Arnold's theorem [T], stating 
that if a one-parameter group of symmetries of a direction field is 
known, the equation dp/dq = f{p, q) can be integrated explicitly. This 
is obvious for the direction field ((71)) of the gauged pendulum. 

3.6 Reduction using Rescriptive Gauge Sym- 
metry 

It is a well-known fact in dynamical system theory that under certain 
conditions, systems that exhibit symmetry are also reducible j3U]. We 
shall discuss reduction in the context of rescriptive gauge theory by 
following a few fundamental steps; ultimately, we will show that re- 
scriptive gauge symmetry allows to reduce classes of nonlinear system 
into linear ODEs, solved by simple quadratures. 

We begin our quest for the manifestation of reduction in the realm 
of rescriptive gauges by asking how a rescriptor for a given ODE can be 
found. We shall then show that the answer to this question is related 
to a more profound problem - that of exact linearization of ODEs, or, 
as we shall call it for clarity - global linearization. We shall dwell upon 
the latter issue shortly, and will first address the more basic query. 

Finding a rescriptive gauge transformation for a given ODE is im- 
portant, since it may allow quadrature in the modified time scale by 
reduction into linear forms. Consider, for illustration, the 1-DOF 
gauged pendulum model 



q = pu{p,q) 
p = -qu{p, q) 



(75a) 
(75b) 
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which is readily transformed into the ODE 



■• dujp, q) 

q — — pq + u{p, q)q 



du{p, q) 
dp 



+ u{p, q) 



(76) 



Thus, any ODE that is written in the form H76|) can be transformed 
into the de-rescribed gauged pendulum q" + q = using the rescriptive 
gauge transformation dr = udt. However, usually the rescriptor, u, 
cannot be easily found. Consider, for instance, the nonlinear ODE 



O 9 

q — q cot q + q sin g = 0, 



(77) 



for which the rescriptive gauge transformation 



dr = sin qdt, 



(78) 



reveals that ()77() is no more than a harmonic oscillator in disguise, 
viz. q" + q = 0. However, one cannot determine that u = sing by 
observation. This calls for a more rigorous methodology for finding 
the rescriptor. 

To that end, consider a second-order ODE of the form 



q + fiq)q'^ + hu{q)q + ijiq) = 0. 



(79) 



When can this ODE be transformed into the linear form 



q" + hq' + boq + c = (80) 

by a rescriptive gauge transformation 

dr = u{q)dt (81) 

only? The answer lies in the theory of exact linearization j^, which 
seeks a transformation rendering a nonlinear ODE amenable for quadra- 
ture. We shall prefer the term global linearization, emphasizing that 
this method is conceptually different from the common point lineariza- 
tion. We shall ultimately use global linearization theory to help us 
track down the rescriptor of a given ODE. 

The theory of global linearization suggests that ODEs of the form 
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()79|) can be globally linearized by a transformation of the from 

z = P J uexp fdq^ dq, dr = u{q)dt 

where (3 = const, if and only if (|8fl|) can be written in the form 

Q + f{q)Q^ + fei-u^ + u exp ^- y f{q)dq^ 
iQ y nexp ^y f{q)dq^ ^ 



(82) 



(83) 



This fundamental result can be adapted to the case in question. In 
particular, since we are probing the case of rescriptive gauge transfor- 
mations, we must require that z = q, or, in other words, that 



P=l, u = u{q), f 



1 du 
u dq 



(84) 



In our discussion we allowed u to be a function of both q and p, while 
()82|) permit a u which is a function of q only. Thus, we must take 
u = u{q), as written in (|84j) . Relations (|84j) modify (jSSl) into 



1 du 2 

— —q + biuq + u exp 
dq 

1 dn 



g —q + 6iug + u exp 

u dq \ J u < 

1 du 2 , 2 7 c 

q + biuq + u 6og + - = 



f{q)dq 
1 dti 



+ 



(85) 
(86) 
(87) 



Thus, we have proven that a second-order ODE can be transformed 
into a linear ODE using a rescriptive gauge transformation (assuming 
that the rescriptor is a function of the coordinate only) if and only if 
this ODE can be written as 



u{q) dq 



Eq. (|88|) immediately yields the rescriptor: It is the square root of the 
coefficient of the coordinate, q, divided by y/bo- 

Returning to example ((77|l , we see that it succumbs to the general 
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from (|S5|) by substituting 



bi=0, 60 = 1, c = 0, (89) 

which yields 

q--^q'^+u^q = (90) 
u aq 

and it is immediately apparent that the rescriptor is u = sing, which 
agrees with H78|) . As a simple verification, we also note that 

1 du 

— - = cot q. (91) 

u aq 

Similarly, the chemical rate equations H59|) . written as the single ODE 
d^[A] 1 fd[A]Y , r,,d[A] 



may be linearized using the transformation dr = [A]dt, as was done 
in E31 

The above process can be repeated for higher-order ODEs as well. 
The bottom line is that the theory of global linearization is a con- 
venient method for finding a rescriptor of a given ODE, or, in other 
words, reduce it into a linear ODE using a rescriptive gauge transfor- 
mation. 

To conclude this section, we shall show that there are well-known 
ODEs that can be transformed into the reducible form (|88|) using an 
additional auxiliary variable transformation. This observation is in- 
spired by where we have shown that the Euler-Poinsot equations 
are transformed into a linear from in the independent variable r us- 
ing the rescriptive gauge transformation dr = u^dt and the auxiliary 
transformation z = a;|. For example, consider the ODE: 

q + qq + kq^ = 0, k = const. (93) 



This ODE arises in a few practical problems [22|. To render it glob- 
ally linearizable using a rescriptive gauge transformation, perform the 
auxiliary variable transformation z = q^,so the modified system reads 
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z - —z^ + ^J^z + kz^ = 0. (94) 

In this form, ()94p adheres to ansatz (|88() . with the rescriptor u = \/z = 
q and k = bo, bi = 1, c = 0. The rescriptive gauge transformation 
dr = y/zdt transforms (|94|) into 

+ z' + 2kz = 0. (95) 



3.7 Illustrative Examples 

We shall now illustrate the rescriptive gauge transformation formalism 
and the resulting gauged pendulum concept using a few numerical 
examples. 

Example 1 (A damped pendulum is a gauged pendulum) 

Consider the model [HS] : 

q = up (96a) 
p = —uq (96b) 
ii = —oJ^q — au. (96c) 

By carrying out the transformation p = cos (f), q = sin (j), these equa- 
tions are immediately recognized as a state-space model for a damped 
nonlinear pendulum, 

<i> + uj'^sm(f) + a4> = 0. (97) 

System (|96]l complies with the gauged pendulum formalism (|26|) : it 
can be therefore viewed as a rescribed harmonic oscillator, revealed 
by the rescriptive gauge transformation d(j) = udt, so that t = (j): 

q'=P, P=-q- (98) 

Obviously, the rescriptor, or gauge velocity, is simply the angular ve- 
locity, i. e. u = (j). The scalar differential equation for this dynamic re- 
scriptive gauge function, Eq. (|96c|) . assumes the nonautonomous form 
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U = —I 



u)q sm[cj))/u — a. 



(99) 



For a = 0, the rescriptive gauge function does not explicitly depend 
upon the rescriptor itself, and (|99|) is easily solved by quadrature: 



It is interesting to note that under the rescriptive gauge symmetry, 
the harmonic oscillator and the damped nonlinear pendulum are rep- 
resented by the same mathematical formalism - although for different 
independent variables - whereas the time flow of these models is com- 
pletely different. The harmonic oscillator, which is a conservative 
system, does not have an attractor, since the motion is periodic. The 
damped pendulum, on the other hand, is a dissipative dynamical sys- 
tems, in which volumes shrink exponentially, so its attractor has 
volume in phase space. This alleged paradox stems from the fact that 
the dissipative time flow of the damped pendulum becomes periodic 
under a change of the independent variable. Thus, an observer mea- 
suring the "time", (/>, is bound to observe periodic behavior, while an 
observer measuring the "true" time, t, will observe exponential decay. 

These observations are demonstrated and validated by means of a 
numerical integration, comparing the flows of ()96() and ()98() . Figure 
n compares between q{t) (Fig. ^) and q{4>) (Fig. ^3); and between 
p{t) (Fig. Ct) and p(</.) (Fig. W), for a = 0.1, go = 0.5, po = 1, ^^o = 
5, (/>o = sin~^ qq = 0.5236. 



2u;q(cos (j) — cos (po) + u'^{(j)o). 



(100) 
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Figure 1: An exponential decay of a damped nonlinear pendulum can be 
transformed into harmonic oscillations by a rescriptive gauge transformation. 

Example 2 (A glimpse of order in the realm of chaos) 

Consider the dynamical system 



q 

P 

X 

y 

z 



yp 
-yq 

a{y - x) 
(r — z)x — y 
xy — bz 



(101a) 
(101b) 
(101c) 
(lOld) 
(lOle) 



where a, r, b are constants. Eqs. HlUlc|) - ()lUle() are recognized as the 
Lorenz system, and the entire system ()1U1() complies with the gauged 
pendulum formalism of Eqs. (|2l-ij) . It shah be thus referred to as the 
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Lorenz-fed gauged pendulum. 

For certain parameter values and initial conditions, the Lorenz 
system is known to exhibit chaos. For instance, choosing the param- 
eter values o" = 10, r = 28, b = 8/3, the initial conditions x{0) = 
10, y(0) = 10, 2:(0) = 10, and simulating foi tf = 50 time units, yields 
the trajectory depicted by Fig. |21 




Figure 2: The Lorenz strange attractor feeding the gauged pendulum. 



Let us now examine the time history of p and q, shown in Fig. |31 
and ask: Do q and p exhibit chaotic behavior? To answer this seem- 
ingly trivial question (without using a comprehensive mapping of the 
phase space using Poincare sections) , we shall resort to the common 
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"engineering" interpretation of chaos, although more mathematically- 
rigorous definitions, related to the destruction of KAM tori |3] or the 
Kolmogorov-Sinai entropy do exist. As Strogatz says in refer- 
ence "no definition of the term chaos is universally accepted yet, 
but almost everyone would agree on the three ingredients used in the 
following working definition" . These three ingredients are: 

1. Aperiodicity: Chaos is aperiodic long-term behavior in a deter- 
ministic system. Aperiodic long-term behavior means that there 
are trajectories which do not settle down to fixed points, periodic 
orbits, or quasiperiodic orbits as t ^ oo. ^ 

2. Sensitive dependence on initial conditions: Nearby trajectories 
separate exponentially fast, i.e., the system has a positive Lya- 
punov characteristic exponent (LCE). 

Strogatz notes that he favors additional constraints on the ape- 
riodic long-term behavior, but leaves open what form they may 
take. He suggests two alternatives to fulfill this: 

3. Requiring that there exists an open set of initial conditions hav- 
ing aperiodic trajectories, or 

4. If one picks a random initial condition x(to) = xq then there 
must be a nonzero chance of the associated trajectory x{t, xq) 
being aperiodic. 

Returning to Fig. El we see that items 001 and |3 in Strogatz's list 
are satisfied: p and q exhibit aperiodic behavior, the open set of initial 
conditions guaranteeing aperiodic trajectories for a = 10, r = 28, b = 
8/3 are xo,yo, zo,po,qo G R\{0}, and hence for randomly selected 
initial conditions p and q will be aperiodic. The only remaining test is 
to calculate the LCEs, denoted by Aj, i = 1, . . . , n. However, as shall 
be illustrated shortly, calculation of the LCEs may be problematic for 
system (fTnT|) . 

First, we should note that some authors endorse the calculation 
of the maximal Lyapunov exponent in order to establish the presence 

^For the purposes of this definition, a trajectory which approaches a limit of oo as 
t oo should be considered to have a fixed point at oo. 
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Figure 3: A seemingly irregular behavior of a gauged pendulum fed by a 
chaotic process. 

of chaos. For example, Ref. |34j states that "it is well-known that 
the ordered or the chaotic property of an orbit is characterized by the 
largest Lyapunov characteristic exponent" . This approach, however, is 
misleading for system (jlOlj) . To illustrate this fact, we have calculated 
the maximal Lyapunov exponent for (jlUlj) using the standard method 
developed by [23] and The result is depicted by Fig. |1] for an 

integration period of 12, 000 time units. It is seen that the maximal 
LCE satisfies maxjAj ~ 0.9, which is the well-known maximal LCE of 
the Lorenz system. Hence, according to the rationale of ,34j, system 
is chaotic - or is it? 
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Figure 4: The maximal Lyapunov characteristic exponent for a Lorenz-fed 
gauged pendulum system. 

For a more rigorous analysis, the entire spectrum of LCEs should 
be examined. Since the LCE spectrum of the Lorenz system is well- 
known (the phase space contraction satisfies the relation ^ - Aj = V • 
[x, y, z] = —((7+6+1) = —13.667), let us concentrate on the additional 
LCEs contributed by p and q. A magnified view of these LCEs is 
shown in Fig. El One of these LCEs is smaller than zero, while the 
other one assumes the value of 4- 10~^, which allegedly indicates that 
the additional states are also chaotic. 



However, this is a mere illusion resulting from the fact that the cal- 
culation process of the LCEs is affected by the truncation and round- 
off errors of the numerical integration routine used to simultaneously 
integrate the extended phase space of the original and linearized sys- 
tems.^ One may view this phenomenon as pseudochaos [221; the truth 



^This causes the Lyapunov exponents themselves to exhibit a chaotic behavior; most 
high-order integrators are chaotic maps, as pointed out in 7 . This may be viewed a 
manifestation of the uncertainty principle. 



31 



X 10" 



-1 

1.191 1.192 1.193 1.194 



1.195 1.196 
Time unit 



1.197 1.198 1.1S 



1.2 



X 10 



Figure 5: The Lyapunov exponents contributed by the states p and q. 

regarding "chaos" in system (jlUip can be plainly revealed by real- 
izing that (jlUip complies with the gauged-pendulum formalism, and 
can hence be subjected to a rescriptive gauge transformation of the 
from dr = ydt. This transformation will transform (|inia|) . ()inib|l 
into q' = p, p' = —q, which is an integrable system and hence cannot 
exhibit chaos. This observation is illustrated in Fig. |H1 showing plots 
of q and p as a function of r. Thus, in contrast to the prediction 
of the common engineering interpreting of chaos and the chaos de- 
tection tools thereof, the rescriptive gauge transformation shows that 
the temporal behavior of signals cannot always be used to predict the 
presence of chaos. This observation calls into being the concept of 
partial chaos 0, meaning that in a given system, both chaotic and 
regular signals may co-exist, even if the chaotic states overshadow the 
regular behavior of the other states. 

Another important conclusion concerns the system observables. 
Observables, or outputs, is a subset of state variables, z, dimz = / < 
dimx = n, determined by the output map, O : M" M', such that 
z = £)(x), and an observation scale, T G R, such that z : R ^ M'. If 
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1. = t, the observation process is temporal and the observable scale is 
merely the time. Our simple example shows that temporal observa- 
tions may be misleading when used to detect chaos, even when using 
a seemingly rigorous test such as the LCE spectrum. A fictitious ob- 
server using T = r as the scale would have not suspected that the 
Lorenz-fed gauged pendulum is a chaotic process. 

We further conclude that rescriptive gauge transformations may 
be used to isolate self-similarities of a dynamical systems. In our 
example, the Lorenz system remains scale-invariant; i. e. its Hausdorff 
dimension does not depend on the scale. However, the Lorenz-fed 
pendulum is not scale invariant, and hence is a regular process in 
disguise. 

Example 3 (Stochastic signals, coding, and Kolmogorov complexity) 

The preceding example illustrated the fact that the gauged pen- 
dulum concept may be used to order pseudochaotic behavior. This is, 
in fact, only an understatement of the potential of rescriptive gauge 
theory; this theory can be used not only for ordering pseudochaotic 
signals, but moreover, transform seemingly stochastic signals into de- 
terministic ones. 

Our final example is therefore concerned with illustrating how re- 
scriptive gauge symmetry, and in particular a simple gauged pendu- 
lum, may be used to establish some key ideas in modern information 
and coding theory through the well-known notion of Kolmogorov com- 
plexity. 

The Kolmogorov complexity (also known as Kolmogorov-Chaitin 
complexity, stochastic complexity, and algorithmic entropy) of an ob- 
ject is a measure of the computational resources needed to specify the 
object |E1 Uni IHl- In other words, the complexity of a string is the 
length of the string's shortest description in some fixed description 
language. It can be shown that the Kolmogorov complexity of any 
string cannot be too much larger than the length of the string itself. 
Strings whose Kolmogorov complexity is small relative to the string's 
size are not considered to be complex. The sensitivity of complexity 
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Figure 6: The seemingly irregular behavior of the Lorenz-fed gauged pen- 
dulum, shown in Fig. [31 can be regularized into harmonic oscillations by a 
rescriptive gauge transformation. 

relative to the choice of description "language" is what the current 
example is about. To that end, consider the gauge pendulum 

q = wp, p = —wq (102) 

where here the rescriptor w is a hand-limited white noise, that is, 
a white noise going through a zero-order hold with some sampling 
frequency Tyj and power spectral density W . Model (|lfl2j) can be 
de-rescribed by dr = wdt. 

Let us compare the representation of the "strings" q and p using the 
"languages" t, time, and r, a random walk obtained by integrating w 
(i. e., a stochastic signal in its own right). This comparison is depicted 
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in Fig.[7|for T^, = 0.1 time units and W = 0.1. Fig. [7^ shows the signal 
q{t), which should be compared to the signal g(r), shown in Fig. [3). 
Similarly, compare p{t), Fig. Et, to p(t), Fig. 




Although q(t) and p{t) seem stochastic and therefore Kolmogorov- 
complex in the "language" t, their alleged complexity vanishes when 
the "language" r is used, and the stormy stochasticity vanishes into 
harmonic oscillations, implying much reduced Kolmogorov complex- 
ity. This phenomenon has practical value in terms of coding theory: 
Signals may be coded using the "code" t and de-coded using the "key" 
r. 
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4 Descriptive Gauge Symmetry 



A more "benign" gauge symmetry can be detected by applying the 
rationale of the Yang-Mills gauge theory of infinite-dimensional sys- 
tems on finite-dimensional dynamics. In this case, the problem can be 
defined as follows: Given a rescribed dynamical system, of the form 
0, 

5't o f = g(x,u(x,t)) = — , 

find a descriptive gauge function ^j, possibly different for each re- 
scriptor component Ui, satisfying 

d^j = Hi(-K,dx,x,dx,Ui(x,t),t,dt), (103) 

such that the gauge automorphism 

5* o V^te(xo, to), t) = V3(Ci(xo, to), t) (104) 

holds for some ^j(x, t), € J^, i S [1, . . . ,n], where is some abstract 
configuration manifold embedded in M". If j £ [l,...,n] sat- 
isfying (|1()3)) such that (|l()4j) holds, then we shall say that system 
(jSJ exhibits descriptive gauge symmetry under the descriptive gauge 
transformation H1U3|) . In this case each - and not the rescriptor, as 
in rescriptive gauge symmetry - becomes either a static or a dynamic 
descriptive gauge function. 

A descriptive gauge symmetry of order k or simply partial descrip- 
tive gauge symmetry comes about when the descriptive gauge trans- 
formation does not affect k state variables, k < n, viz. 

5* o(^(^i(x,to),t) = ^(ei(x,to),t), (105) 

4.1 Newtonian Systems Revisited 

Consider the Newtonian system 

q(t)+f[q(t)] =u[q(t),q(t),t] (106) 
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with the configuration manifold Q, q G Q, and the tangent bundle 



TQ 
Q X 



Bn/2 



Bn/2 



, SO that (q,q) G TQ, f : Q ^ M"/^ and u 



is the rescriptor. Let 



q = 7[x(t),t] 



(107) 



be the solution of (|106|) . where x : M ^ C M" are the variational 
coordinates. The velocity vector field is then given by the Lagrangian 
derivative® 



q = /3[x(t),x(t),t] = ^^[-5'^U^^[y^'^]x, 



dt 



(?X 



(108) 



or, stated in terms of field theory, the gauge covariant derivative of 
the configuration vector field. 



Dt7 = da + Vx7 • Dtx, 



(109) 



and the velocity vector field, 



Bt(3 = dt(3 + Vx/3 • Dix. 



(110) 



Denote the convective (sometimes also called advective) term by 



dx. 



(Ill) 



Substituting (fW]) and (|TnH|) into (|Tn?T|) transforms (PT7H|) into the 
Gauss-Poisson equations 



dt 



P^(x) 



^7 



-I T 



d-^\ 



dx. 



-I T 



(112) 



where P is the n x n skew-symmetric Poisson matrix, whose entries 
are the Poisson brackets, {xi,Xj}. 

We note that q - the physical trajectory on the configuration man- 
ifold - remains invariant under any selection of ^. Thus, we are in 
the liberty of choosing a descriptive gauge function vector of the form 



derivative taken with respect to a moving coordinate system. Alternatively, this 
operation is sometime referred to as the substantive derivative or Stokes derivative. Fluid 
dynamicists prefer the notation d/dt = d/dt + v • V, where v is the velocity vector field. 
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^ ^ I W(x, u(x, t),±, t), u(x, t) + 
1 0, u(x,t) =0, 

inducing an affine descriptive gauge transformation of the form 

= H(x, dx,x, dx, (114) 

with 

F(x, dx, X, dx, t) = dt^ + Vx* • X + Vi* • x, (115) 
such that the gauge automorphism 

o 93(q(xo, to), i) = V'(q(xo, to), t) (116) 

holds. We note that the gauge ^ in the Newtonian context has dimen- 
sions of velocity, and hence can be referred to as the gauge velocity. 
Recall that we have made a similar observation regarding rescriptive 
gauge symmetry in Newtonian systems (cf. ^3.21) . 

We see that Newtonian systems exhibit partial descriptive gauge 
symmetry, so that trajectories in the configuration space remain in- 
variant under a selection of a particular descriptive gauge function. 
Stated more eloquently, q = 7[x(t),t] remains invariant under the 
symmetry transformation 

da ^ Dn = da + (ii7) 

The gauge group Q therefore consists of real valued functions on R"/^ , 
with the group operation being addition. An element ^ acts on the 
velocity vector field according to the rule H117() . 

Eq. ()112|) is not necessarily integrable, and may posses no "classi- 
cal" integrals whatsoever. However, regardless of the particular prop- 
erties of the original system (|106)) . the variational system (|112)) must 
always satisfy the constraint H118|) . If we choose = 0, then H118|l 
becomes 

d-f\x(t),t] 

^^r} ' ^ x = 0, 118 
ax 

which may be viewed as a hidden integral emanating from the descrip- 
tive gauge symmetry. 
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We emphasize that our hidden symmetry is not confined to systems 
in which the homogenous solution of q + f(q) = can be found. 
Although there are important realms of science in which a solution 
to this system does exist - the most notable being the case where 
f = Vi?, where R is an inverse square gravitational potential emerging 
in Keplerian orbital mechanics - in many other instances -f cannot 
be found in closed form. This stems from the fact that the distinction 
between f and u, the rescriptor, is really an artificial one; we can 
always take % = u — f , so that now 

q(t) =x[q(t),qW,t]. (ii9) 

Letting x = [x^ entails 

7[x(t),t] =xii + X2. (120) 
In this case the Gauss-Poisson equations are simply 



Xl 




0n/2xn/2 ^ 






X2 




I -tl 




X(xi,X2,*,t)-* 



(121) 



where / is an n/2 x n/2 identity matrix and 0^/2 xn/2 is an n/2 x n/2 
zero matrix. 

For ^' = 0, our hidden integral re-emerges, assuming the particu- 
lary simple form 

xit + X2 = 0, (122) 

for which 

Xl = x(xi,X2,t) (123a) 
±2 = -tx(xi,X2,t). (123b) 

Eqs. (|123|) are particularly amenable for numerical integration, be- 
cause in this form constraint (|122|) should be satisfied; however, due 
to numerical round off errors, this constraint is violated. An improved 
numerical integration may be achieved if the integration scheme itself 
is forced to satisfy this constraint during the integration process. 
We shall illustrate these observations in ^4.31 discussing a few nu- 
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merical examples showing how descriptive gauge symmetry may be 
used to improve numerical integration of ordinary differential equa- 
tions. Our last section before dwelling upon actual examples deals 
with reduction in descriptive gauge theory. 

4.2 Reduction using Descriptive Gauge Sym- 
metry 

Equivalently to ^3.61 we shall conceive a process for reduction using 
descriptive gauge symmetry. To that end, we re-write Eq. (|112p into 
the following form: 

n/2 n/2 

Xj{t) = '^fji{x){ui-ipi) -'^gji{x)Tpi, j = l...n (124) 

i=l i=l 

This yields n/2 integrals x^, k G N"/^ (there are n!/(n/2!)^) pos- 
sible combinations of constants of motion obtained by concomitant 
n!/(n/2!)^ descriptive gauge function components) obtained by solv- 
ing the n/2 first-order ODEs 

n/2 n/2 n/2 

i=l i=l i=l 

The freedom to reduce system (|112|) stems from the existence of the 
descriptive gauge function ^. If we fix the gauge - a straightforward 
selection would be ^' = 0, as Lagrange himself had advocated in 
his memoirs [211 1221 I23j - this freedom will be lost, and hence the 
possibility for reduction. This process is illustrated in the following 
section. 
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4.3 Illustrative Examples 

Example 4 (Reduction using descriptive gauge symmetry) 

Our first example is a simple one, illustrating the concept of re- 
duction using descriptive gauge symmetry. To that end, consider the 
one-dimensional, second-order ODE 

q{t) + q{t)=sm{t), (?(to) = 9o, q{to) = %■ (126) 

The general solution is 

q{t) = xi{t)qi{t) + X2{t)q2{t), (127) 

where qi{t) and q2{t) are the fundamental solutions 

qi{t) = cos(t), q2{t) = sin(t). (128) 

Taking ^' = ^'(t) to be a time-dependant descriptive gauge function, 
the Gauss-Poisson equations (|112)) assume the simple form 

xiit) = ^it)q2it) - q2it)[uit) - 4>{t)] (129a) 



X2{t) = qi{t)[u{t) - ^(t)] - ^{t)qi{t) (129b) 

Substituting (|128|) into (|129|) yields ODEs for two possible descriptive 
gauge functions that will transform system H129() into a single ODE. 
Thus, letting ki and /c2 denote arbitrary integration constants, taking 
$(t) = $i(t) = [t/2 - sin(2t)/4 + ki]/ sint will reduce (TT^ into 

xi = const. (130a) 
±2 = -l/2[costsint - t - 2A;i]/sin^t (130b) 

and taking ^'(t) = ^'2(t) = [-cos(2t)/4 + fc2]/cost wiU reduce ([T^ 
into 

X2 = const. (131a) 
xi = -l/4[2cos'^t-l-4k2)/cos'^t. (131b) 
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Both H130b|) and H131b|) are readily solved by quadrature. Systems 
(|13()|) and (ll.Slj) will of course both yield the same general solution 
(|127|) : this is what descriptive gauge symmetry is all about. 

Example 5 (Gauss-Poisson variables reduce Hamiltonian drift) 

Consider the Hamiltonian 

n=l{q' + q') (132) 

of the harmonic oscillator 

q{t) + q{t) = 0. (133) 

We shall compare the numerical integration of this equation in two 
cases. In the first case, the state variables are chosen in standard 
form: qi = q, q2 = q, so that p33jl becomes 

qi = q2, q2 = -qi, (134) 

while in the second case, ()133() is re-written in the form q = —q, and 
the solution is taken as 7 = xi{t)t + X2{t), with gauge ^' = 0. Per 
(|123|) . this yields the state-space representation 

Xi = —Xit — X2, X2 = Xit^ + X2t. (135) 

We used MATLAB's ODE45 integration routine, a 5th-order Runge- 
Kutta integrator with an adaptive time step, to integrate H134() and 
()135() with an integration tolerance of 10~^ for 5000 time units given 

n = 2.5. 

The time history of the Hamiltonian for both cases is plotted in 
Fig. ISl Since the ODE45 routine is not a symplectic integrator, the 
standard selection of states, Eq. ()134|) . causes the Hamiltonian to de- 
crease with time at a rate of about 5 • 10""^ units per time unit. This 
introduces artificial numerical damping into the system, which causes 
the harmonic oscillations to slowly damp out. However, the same sys- 
tem integrated in form H135() keeps the Hamiltonian fixed. This implies 
that the Gauss-Poisson formalism may be used to "symplectify" non- 
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symplectic integrators by a judicious selection of the descriptive gauge 
function. 
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Figure 8: Symplectifying a non-symplectic integrator using Gauss-Poisson 
state variables and a zero descriptive gauge function. 



Example 6 (Descriptive gauge reduces integration errors) 

Our final example shows how descriptive gauge symmetry may me 
used to reduce the numerical truncation error of ODE integration. 
We shall ultimately show that using the Gauss-Poisson state variables 
with an appropriate descriptive gauge can dramatically reduce the 
numerical truncation errors, and show how such gauge can be found. 
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Consider, for example, the one-dimensional, second-order ODE 



qit) + 2^uJnqit) + ujlqit) = n(t), q{to 



), q{to) 



(136) 



where ujn is the natural frequency, is the damping coefficient and the 
rescriptor u{t) is piecewise continuous. Assuming an underdamped 
case < 1), 



q{t) = Xi{t)qi{t)+X2it)q2it), 
where qi{t) and q2{t) are the fundamental solutions 

= e-«""*cos(u;rfi) 



(137) 



(138a) 



92 (t) 



sm{uJdt) 



(138b) 



and u>d = u>n\/'^ — i"^- Taking ^! = ^{t) to be a time-dependant 
descriptive gauge function, the Gauss-Poisson equations (|112j) assume 
the simple form 



xi{t) 



^{t)q2{t) - q2{t) U{t) - i!{t) - 2iu:n^{t) 



w[qi{t),q2{t)] 



(139a) 



X2{t) 



qi{t) u{t)-i!{t)-2iujn^{t) -^{t)qi{t) 



where w[qi{t) , q2{t)] is the Wronskian determinant. 



qiit) q2{t) 
qiit) q2{t) 



The initial conditions for this system are then 

-qm2{to) + q2{to)qo - 92(^0)^(^0) 



xi{to) 
3^2 (to) 



91(^0)92(^0) - 92(^0)91(^0) 
-9091(^0) + 91(^0)90 - 91(^0)^(^0) 
91(^0)92(^0) - 92(^0)91(^0) 



(139b) 



(140) 



(141) 
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From the discussion in ^4.11 we know that there is only a single solu- 
tion for q{t) for any given initial conditions; thus, per the descriptive 
gauge symmetry, q{t) must remain invariant to any selection of the 
gauge function ^. However, ^ may be used as a tuning function for 
mitigating the numerical integration error. To that end, we define the 
numerical integration error of some state variable (•) as the difference 
between the true solution and the numerical solution: 



{')true {,') numerical 

(142) 

In this example we will demonstrate how to mitigate the numerical 
integration error of a (fixed-step) 4th order Runge-Kutta integrator 
(RK4) by several orders of magnitudes. This merit is achievable by 
applying gauge- optimized integration. To that end, let us re-write 
Eqs. (fT39|l into 

xi{t) = f{t,^) (143a) 
X2it) = g{t,^) (143b) 

Obviously, in the linear case discussed herein, the transformation 
into the Gauss- Poisson equations has transomed the ODE integration 
problem into a simple quadrature, whose accuracy can be controlled 
by a proper selection of a time-dependant descriptive gauge function. 

The integration errors resulting from numerically integrating (|143|1 
are given by [33 

e., = -^h'f^'^[^,^m (144a) 
e., = -1/^5^(4) [e,^(e)] (144b) 

where G (Iq, t). The total integration error of q{t) is now calculated 
as follows: 



eq = qi{t)ex^ + q2{t)ex2 

= Qiit)[-^h'f^'\^^n+Q2{t)[-^h'g'^'\^,n- (145) 
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If some 4** could be found for which eg = 0, G {tQ,t), then the only 
remaining integration error of q{t) would be the numerical round-off 
error. Indeed, such can be quite straightforwardly found. For 
example, if the forcing term is of the form 

F{t) = h Ofe cos{kujot) + bk sm{kuJot), (146) 

fc=l k=l 

then ^* may be chosen as a Fourier series as well: 

*'(*) = ^ + ^[An cos{nujot) ]+Y.[^^ sin(na;oi) ] (147) 

n=l n=l 

Substituting (|147|) into (|145|) and solving for the coefficients by requir- 
ing Cq = yields, after some algebra, that 

= . X i^n [an cos{ntujo) - bn sin(ntc<;o)] 

+ ——-ujni [a-n sin(nta;o) + bn cos{ntuQ)\ (148) 
^-^ pin) 

n=l ^ 

where 

pi = 4[32u;U^ + (-40^^ + Sln^ivt^l)^^ + {-1^r?^i^l + 14a;^ + lOn^u;^^^)^' 

P2 = -4n(16(^6^6^(16n2c^^u;2-16u;n)e' + (-12n2^^u;o' + 4a;n)e' 
+5u;gn6 + lln^c^^wg + ^6 ^ ISn^w^o;^) 

p = {2mu)lf + (320a;^wgn2 - 384u;^)^'^ + (176a;^ 160u;^a;^n^ - 'i2Qujli^lr?)i'^ 
+100a;^w^n^-h25a;gn^). 

(149) 

For illustration, if we desire to integrate numerically 

q + q = s\ii2t, g(0)=0, g(0) = 0, (150) 
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then the descriptive gauge function yielding minimum numerical trun- 
cation will be, based on (|14<S|) . simply 

^'^ = -^cos(2t). (151) 

In Figure IHl we depict a comparison of integration errors between 
the gauge-optimized integration utilizing the Gauss- Poisson equations 
with the optimal descriptive gauge function (jlSlj) and the standard 
choice of state variables qi = q,q2 = Q- As can be plainly seen, the 
gauge-optimized integration decreases the integration error by three 
orders of magnitude in the examined time interval. Moreover, the inte- 
gration error using the standard state variables is diverging, while the 
error of the gauge-optimized integration is bounded. Therefore, for 
a larger time interval, the use of gauge-optimized integration, utiliz- 
ing the concept of descriptive gauge symmetry, becomes increasingly 
important. 



5 Summary and Conclusions 

This paper described how gauge theory can be adapted for finite- 
dimensional dynamical systems. We have defined gauge symmetry 
in a very broad context, and distinguished between two fundamental 
manifestations of gauge symmetry: 

(i) Rescriptive gauge symmetry results from an action of a one- 
parameter Lie group, yielding an Abelian Lie algebra. A rescriptive 
gauge symmetry transformation is then an infinitesimal change of the 
independent variable, which renders the system integrable via reduc- 
tion. 

(ii) Descriptive gauge symmetry is an invariance of some configu- 
ration space under a gauge transformation of the covariant derivative. 
In this case the symmetry group consists of real- valued functions on 
the Euclidean space, with the group operation being addition. 

The gauge conversation leads to a few practical conclusions. We 
first note that gauge symmetry is ubiquitous in a myriad of scientific 
fields. Gauge theory for finite-dimensional system may be thus viewed 
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Figure 9: The numerical integration error of the gauge-optimized integration 
considerably reduces the integration error compared to a standard choice of 
state variables. 

as a generalization of dynamical systems theory into the realm of group 
theory, unifying various physical phenomenon into simple generating 
models. 

Furthermore, the gauge-theoretic tools may be used to improve our 
understanding of chaos, randomness and their inter-relations. We dis- 
cussed a few simple examples showing how a change of scale can lead 
to pattern evocation in seemingly chaotic and/or stochastic systems. 

Finally, gauge-theoretic tools are important for improving the ac- 
curacy of numerical integration. The gauge freedom allows re-shaping 
of the phase space so as to render it tractable for numerical integra- 
tion. 
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